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ABSTRACT 


The  aerodynamic  roughness  of  the  extremely  complex 
terrain  at  White  Sands  fissile  Range  has  been  examined 
through  the  analyses  of  wind  and  temperature  profiles 
measured  on  the  62 -meter  meteorological  research  tower. 

The  roughness  length  zQ  was  considered  as  a  function  of 
the  season  and  the  direction  of  the  mean  flow.  There 
was  an  arparem:  correlation  between  the  roughness  length 
and  the  growing  season  with  variations  in  zq  from  approxi¬ 
mately  10  cm  in  January  to  22.5  cm  in  July.  Variations 
of  nearly  60  cm  were  noted  when  z  was  considered  on  a 
directional  basis. 
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INTRODUCTION 


The  determination  of  the  aerodynamic  roughness  of  the 
surface  in  the  immediate  vicinity  of  a  data  collection  site 
is  a  basic  requirement  cf  any  study  pertaining  to  the  tur¬ 
bulent  characteristics  of  the  planetary  boundary  layer. 

An  aerodynamically  rough  surface  may  be  defined  as  one  on 
which  the  irregularities  are  large  enough  to  prevent  the 
formation  of  a  laminar  sublayer  so  that  the  mean  flow  is 
effectively  turbulent  down  to  the  surface  (Schlichting, 

I960).  If  the  surface  is  homogeneous  with  a  random  dis¬ 
tribution  of  roughness  elements,  the  roughness  length  may 
be  determined  from  surface  drag  measurements,  observed 
neutral  profiles,  or  thermally  stratified  data  and  appli¬ 
cation  of  diabatic  surface  boundary  layer  hypotheses.  The 
roughness  length  is  an  indicator  of  the  effective  aerody¬ 
namic  roughness  of  a  surface  defined  as  the  height  where 
the  velocity  profile  extrapolates  to  zero.  Unfortunately, 
a  homogeneous  surface  that  allows  the  development  of  equi¬ 
librium  flow  from  all  directions  and  during  all  seasons 
is  rare;  thus  the  aerodynamic  roughness  becomes  a  function 
of  the  local  terrain  and  the  vegetive  canopy,  making  the 
roughness  length  difficult  to  determine  on  a  representative 
basis.  Experimental  data  must  be  carefully  analyzed  in 
terms  of  small-scale  local  terrain  features  and  in  terms 
of  the  growing  season. 

The  purpose  of  this  paper  is  to  examine  the  roughness 
length  of  the  terrain  in  the  .cinity  of  the  62-meter  met¬ 
eorological  research  tower  at  White  Sands  Missile  Range  as 
a  function  of  season,  wind  direction,  and  the  vegetive  can¬ 
opy,  The  terrain  in  the  tower  locale  is  extremely  complex 
and  not  homogeneous  in  all  directions,  but  is  considered  typical 
for  WSMR,  Consequently,  the  mean  flow  is  net  in  equilibrium 
with  the  surface  at  all  times, 

SITE  DESCRIPTION 


The  tower  is  located  in  an  area  that  consists  of  small 
brush-covered  hillocks  one  to  three  meters  hi >*h  and  randomly 
distributed.  Patches  of  sand,  gravel,  native  vegetation, 
i.e.,  bunch  grass,  pigweed,  and  Russian  thistle  (tumbleweeds) 
are  interspersed  between  hillocks.  The  brushy  portion  of  the 
vegetation  consists  of  mesquite,  greasevood,  and  desert  willow. 


Modifications  to  the  natural  terrain  such  as  clearing, 
leveling,  construction  of  buildings,  etc.  contribute  to  the 
overall  heterogeneity  of  the  surroundings  since  the  tower 
is  located  in  an  area  of  extensive  human  activity,  ft.  por¬ 
tion  of  the  natural  terrain  west  of  the  towei  remained  un¬ 
disturbed  and  provided  a  quasi-homogeneous  fetch  of  approxi¬ 
mately  It  kilometers  from  220  degrees  through  350  degrees, 
during  the  first  data  collection  period  from  April  1956 
through  April  I960.  The  determination  of  the  roughness 
length  has  been  restricted  to  the  times  when  the  mean  wind 
was  across  this  130-degree  sector. 

Modification  to  the  terrain  in  this  sector  occurred 
during  the  second  data  collection  period  from  January  1962 
through  October  1963.  Approximately  two  acres  were  cleared 
and  leveled  west  of  the  tower  site,  giving  an  85-meter  fetch 
across  bare  sand  to  the  natural  terrain  and  vegetation.  This 
new  roughness  discontinuity  altered  the  mean  flow  turbulent 
characteristics  of  the  experimental  site  drastically. 

Data  collected  in  the  first  observational  period  are 
representative  of  equilibrium  '‘low  over  a  quasi-homogeneous 
surface,  whereas  those  of  the  second  period  possess  non¬ 
homo  geneous  turbulent  characteristics. 


THE  ROUGHNESS  LENC^H 


Sutton  (1953)  states  that  the  roughness  length  of  a 
surface  may  be  approximated  from 
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(1) 


where  z0  is  +he  roughness  length,  z  the  average  height  of 
the  roughness  elements,  and  C  is  a  constant,  usually  taken 
as  numerica1lv  equal  to  }0  based  on  the  work  of  N'ikuradse 
(See  Sutton,  1953,  and  Schl ichting,  I960). 

For  conditions  of  'V.Hy  rough  turbulent  flow  in  the 
surface  boundary  layer  o:  the  atmosphere  in  neutral  strati¬ 
fication,  the  differential  equation  for  the  velocity  profile 
may  be  writ ‘en  as 
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where  V  is  the  mean  horizontal  velocity,  u*  the  dynamic 
friction  velocity,  k  is  Karman’s  constant  and  z  is  height. 
Integration  yields 
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u* 
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where  z  becomes  a  constant  of  integration.  Over  extremely 
rough  surfaces,  those  covered  with  high  vegetation,  i.e., 
cereal  -'rains  or  trees,  the  velocity  profile  becomes  mere 
complex  and  (3)  does  not  held.  An  empirical  correction  in 
the  form  cf  a  zero-nlane  disDlacement  factor,  d,  can  be 
introduced  in  the  form 
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Equation  (li)  may  be  stated  as  (Lemon,  196? ) 
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where  z  is  the  geometric  near,  cf  the  levels  cf  interest  and 

a 

D  *  d  +°z  ,  the  "effectiveness  displacement  parameter". 
Equations  (3),  (M>  and  (5)  ray  be  solved  graphically  or 
algebraically  for  numerical  values  of  D,  d,  and  7Q  if  repre¬ 
sentative  neutral  wind  profiles  are  availalle. 

I'  the  etrosrhero  is  * herrally  stratified,  equ.at  icr.s  (■) 
*  trough  (3)  no  longer  .  rovide  a  valid  sol'*.',  on  for  ♦  he  wind 
profile.  It  has  been  demon  s*rn*>  d  :y  l'o~;n  a:;-:  .-.krc-v  (If-.!, 
^urinver  (19^1',  rano f r'<?  C 1  °f •  3 '  ar.d  c.l.-rs,  *'.->*  ;•  d.a:a*io 
renditions  the  wind  cir  lx  ronr»>:  or*  od  tv 
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where  ^(z/L’ )  is  a  universal  function  and  L*  is  the  Panofsky, 
ELackadar,  and  KcVehil  (i960)  gradient  length.  Equation  (6) 
predicts  a  straight  line  if  the  velocity  is  plotted  as  a  func¬ 
tion  of  In  z  -  with  zQ  as  the  intercept. 

The  universal  function  $(z/L* )  is  a  stability  correction 
factor  that  has  been  developed  empirically  outside  the  basic 
structure  of  the  Monin-Obukhov  Similarity  Theory  to  compensate 
for  thermal  stratification  effects  on  the  wind  profile.  A 
reasonably  good  fit  to  experimental  data  can  be  obtained  with 
Equation  (6)  in  the  stability  range  0.08  >  Ri  >  -0.10  (at  a 
reference  height  of  1.5  meters)  where  Ri  is  the  gradient  form 
of  the  Richardson  number. 

Restrictions  on  the  use  of  Equations  (3),  (5),  and  (6) 
for  graphical  or  algebraic  determination  of  the  roughness 
length  for  an  experimental  data  site  depend  on  the  thermal 
stratification  of  the  boundary  layer,  surface  homogeneity, 
or  lack  of  it,  stationarity  of  the  mean  flow  and  growth  or 
decline  of  the  vegetive  canopy.  With  these  limitations  and 
restrictions  in  mind,  experimental  data  measured  on  the  62- 
meter  research  tower  were  used  to  determine  the  roughness 
length  of  the  typical  terrain  of  White  Sands  Missile  Range. 


THE  AERODYNAMIC  ROUGHNESS  OF 
AN  EXTREMELY  COMPLEX  TERRAIN 


During  the  25-month  period  from  April  1958  to  April 
I960,  l6ll  profiles  of  wind  and  temperature  for  the  first 
62  meters  of  the  planetary  boundary  layer  were  collected. 

The  raw  data  consisted  of  10-minute  averages  from  nine  levels 
on  the  tower  observed  eight  times  a  day  (0100,  OljOO,  0700, 
1000,  1300,  1600,  1900,  and  2200  hours  MST).  Analyse  of 
these  data  by  Carnes  (1961),  Tourin  and  Hoidale  (1961),  and 
Swanson  and  Cramer  (1965)  indicated  that  a  majority  of  the 
samples  had  been  observed  under  conditions  that  were  suffi¬ 
ciently  stationary  and  homogenepus  to  permit  detailed  analy¬ 
sis  of  the  data  to  establish  the  aerodynamic  roughness  of 
the  surface. 

Owing  to  the  nature  of  the  terrain  and  the  effect  of 
vegetation  on  the  apparent  surface  roughness,  the  analysis 
considered  the  season  and  the  mean  wind  direction  across 
the  130-degree  sector  to  compensate  for  the  homogeneous 
surface.  Profiles  observed  during  the  afternoon  hours  were 
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treated  with  caution  owing  to  the  possibility  of  induced 
error  (KcCready,  1966). 

Usable  profiles  were  first  classified  by  stability 
regime,  utilizing  the  nine  regime  scheme  of  Hansen  (1966). 
Richardson  numbers  for  each  profile  were  computed  from 
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where  g  is  the  acceleration  of  gravity,  T  the  mean  tempera¬ 
ture  in  degrees  Kelvin,  u  and  v  the  north-south  and  east- 
west  components  the  mean  wind,  T  the  adiabatic  lapse  rate 
(“0.0098  deg  km”1),  and  z  the  geometric  mean  of  the  layer 
z  -  z  where  the  subscripts  indicate  the  level  of  interest. 
Only  triose  profiles  where  the  value  of  Ri  was  within  the 
limits  0.08  >  Ri  >  -C.li  at  7.U  meters  were  used  in  the  analy¬ 
sis,  with  emphasis  on  neutral  profiles  when  possible.  Graph¬ 
ical  solutions  to  Eouations  (3),  (5),  and  (6)  were  fitted  to 
the  experimental  tower  data  for  observations  below  approxi¬ 
mately  16  meters.  Above  this  level,  direction  can  no  longer 
be  considered  invariant  with  height,  and  the  approach  to  the 
geostrophic  wind  enters  into  the  solution.  Stability  cor¬ 
rections  to  the  diabatic  profiles  were  obtained  from  Hansen 
(1966a)  and  are  based  on  the  KEYPS  function  which  is  consider¬ 
ed  the  most  reliable  and  simplest  to  use  of  the  Similarity 
Theory  log-linear  profile  hypotheses. 


Average  values  of  the  roughness  length  for  the  research 
tower  site  by  time  of  year  and  wind  direction  are  shown  in 
Figures  1  and  2,  There  is  a  correlation'  between  seasonal 
values  and  the  growing  season,  with  minimum  z  occurring  in 
January  and  a  maximum  in  July  which  represents  the  rainy  sea¬ 
son  for  the  area.  The  directional  dependence  of  the  rough¬ 
ness  length  correlates  very  well  with  the  isoline  analysis 
of  the  turbulence  intensity  measurements  of  Horn  and  "ravle 

(196I4). 


The  second  observational  period  of  1962-1963  gave  ap¬ 
proximately  2000  hours  of  data  that  are  considered  to  be  non- 
homogeneous  owing  to  the  roughness  discontinuity  created  by 
the  sandy  area  west  of  the  tower  site.  Sutton  (1953),  Priestley 
(1959),  and  Lumley  and  Panofsky  (196U)  show  that  the  roughness 
of  bare  sand  should  yield  a  zQ  value  on  the  order  of  0.1  cm, 
while  Sellers  (1965)  indicates  z0's  in  the  neighborhood  of 
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0,03  to  0.0?  cm  for  a  similar  soil  surface  near  Tucson, 
Arizona.  Application  of  graphical  profile  analysis  tech¬ 
niques  to  selected  samples  of  the  .1962-1953  data  gave  a 
roughness  length  of  l."?9  cm  with  a  ±0.79  cm  RKS  error 
which  is  indicative  of  the  heterogeneity  of  the  mean  flow 
across  the  roughness  discontinuity.  Using  the  approach 
of  Dhiiip  (1959),  I>  ’nr  (196}),  Rider,  -'hilip,  and  Pradlev 
(1963),  and  Panofsky  and  '’’owns  end  (196L)  it  was  found  that 
at  1.5  meters  above  ‘he  surface,  the  vine  profile  was  less 
than  two  percent  adjusted  to  the  new  surface  after  travers¬ 
ing  the  C5-meter  fe^ch.  It  is  clear  that  the  nonequilicrinm 
conditions  existing  at  the  nlane  of  measurement  resulted  in 


the  erroneous  1.99  cm  roughness  length  for  the  bare  sand 


fetch  as  compared  with  the  inferred  value  of  0,1  cm. 


DISCUSSION 


Panofsky  (1963)  has  c.emors' rated  the  effect  of  disturbing 
the  natural  vegetation  of  an  area  ,,'^en  the  experimental  site 
was  mowed  during  the  Great  Plains  Turbulence  Study  (See  Lettau 
and  Davidson,  1957),  changing  the  roughness  length  approxi¬ 
mately  30  percent.  Halstead  (1957)  has  shewn  that  the  surface 
roughness  at  the  O'Neill,  Nebraska,  site  of  Project  Prairie 
Grass  (See  Haugen,  1958,  and  Barad,  1 959)  varied  from  0.5  to 
1.1  cm  owing  to  changes  in  vegeo  tion  avid  to  differences  in 
wind  direction  across  a  cuasi-homogeneous  terrain.  These 
causes  and  effects  exist  at  the  WSNit  tower  site.  The  change 
of  surface  roughness,  as  a  function  of  wind  direction  and 
season,  is  directly  attributable  +c  the  grow*1'  and  defoliation 
of  the  vegetive  canopy  and  a  nonuniforr.  terrain. 

The  current  investigation  reveals'  that  the  terrain  in 
the  vicinity  of  the  tower  is  extremely  rough  and  not  homogene¬ 
ous  with  wind  direction  or  seasons,  ^he  mean  of  all  rough¬ 
ness  lengths  determined  fer  the  site  was  17. Ii  cm.  An  inspec¬ 
tion  of  Figures  1  and  2  reveals  that  this  is  a  reasonable 
value  for  the  annual  roughness  regime,  but  could  be.  unreal¬ 
istic,  especially  when  the  directional  data  are  considered, 

The  best  fits  +0  experimental  data  are  obtained  by  use  0^  the 
directional  roughness  length  values,  a] though  the  annual 
regime  curve  will  yield  acceptable  results  within  the  limit 
of  experimental  error  (+  10  per  cent).  For  wind  directions 
other  than  those  within  the  130-degrec  sector,  it  was  fr'nd 
that  a  reasonable  estimate  for  z  was  approximately  20  cm 
which  rives  solutions  within  1  28  rer  cent. 


o 


7*  was  noted  that  a  zero  plane  oisplacem.ent  correction 
to  the  experimental  data  -.-.•as  not  necessary  ever,  during  the 
periods  in  which  The  brushy  portion  of  the  vegetation  was  in 
"ull  leaf.  It  arrears  that  the  hardy  desert  bush  is  flexible 
to  the  roint  that  little  or  nc  bendinr  occurs  even  with  brisk 
wind  conditions. 


The  terrain  notification,  which  occurred  before  the 
.1962-1963  data  collection  period,  resulted  in  a  roughness 
di scontinuity  of  3 rent  100-220  to  1,  approximately  £$  meters 


west  0: 


he  t  ov.'er. 


Ar 


Lnvesti ~ation  of  ’wind  and  temperature 


profiles  observed  ’.’her.  rhe  rear,  '’low  vres  across  the  discon- 

/-• - ---  •• - pofr)  s’novec>  *hat  "kinking"  of 

one  a  cistinct  internal  r.cmentur. 


vinu-ty  \  .p.r.ser.  anc  .-.an: 
"the  c^curr* 


boundary  form.ee  with  an  interface  arrrcxinately 
above  the  surface.  mh:s  agrees  with  the  hypotheses  of 
Elliott  (1956)  and  fanofslc/  and  To’-aisend  (l?6u)  that  predict 
the  interface  to  have  a  height-to-fetch  ratio  of  1:10. 


The  abrupt  roughness  c nance  has  little  effect  on  the 
magnitude  of  the  lateral  and  longitudinal  intensities  of 
turbulence  observed  at  the  tower  as  noted  by  denser  (l?66) 
who  found  few  differences  in  these  parameters  when  they  were 
compared  with  values  of  Swanson  and  Cramer  (1969)  for  the 
195k  «1?60  data.  The  major  contribution  of  surface  discon¬ 
tinuity  was  the  generation  of  nonequi librium  flow  and  ths 
internal  momentum;  boundary. 


CONCLUSIONS 

The  aerodynamic  roughness  of  an  experimental  site  may 
be  determined  from  the  analysis  of  wind  and  temperature  pro¬ 
files  if  the  mean  flow  is  considered  to  be  stationary  and 
homogeneous  across  a  sufficiently  long  uniform  fetch  to  in¬ 
sure  equilibrium,  conditions.  The  growth  of  vegetation  and 
slight  nonuni f ormi ties  in  the  surface  with  respect  to  wind 
direction  are  other  considerations.  The  surface  roughness 
of  a  sparsely  vegetated  and  extremely  uneven  terrain  can  be 
determined  if  caution  is  exercised  ir  the  analysis  of  +he 
profile  data. 

Modifications  to  the  natural  terrain  nave  beer  shown 
to  arfect  the  aerodynamic  characteristics  of  a  surface  to 
the  extent  t.ln+  the  mean  flow  is  no  lenrer  in  eouili braur* 
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wi+h  the  underlying  surface.  The  res' ’’tant  nonhcnsogeneous 
turbulent  flow  completely  negates  the  use  of  conventional 
profile  analysis  techniques  for  determining  surface  rough¬ 
ness  and  forces  the  use  of  estimated  values.  The  results 
presented  m  viyures  1  and  2  are  considered  to  be  represent¬ 
ative  of  the  roughness  length  in  the  vicinity  of  the  62- 
meter  research  tower,  whereas  the  estimate  of  O.I  cm  for 
+  he  65 -muter  bare  fetch  can  be  in  error. 


‘H19N31  SS3NH900& 


FIGURE  2.  THE  AVERAGE  ROUGHNESS  LENGTH  BY  DIRECTION  FOR  THE  SUMMER,  SPRING,  FALL  AND 
WINTER  SEASONS  FOR  WHITE  SANDS  MISSILE  RANGE,  NEW  MEXICO. 
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